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[10] S. Werner and J. Lilleberg, "Downlink channel decorrelation in CDMA 
systems with long codes," in Proc. of 49th VTC, pp. 1614-1617, 1999. 

[0002] This invention pertains in general to communication systems. 
More particularly, embodiments of the invention pertain transmit diversity and 
5 Multiple-ln, Multiple-Out (MIMO) transmission and receiving methods for 
multiple antenna technology of Code Division Multiple Access (CDMA) type 
systems. 

[0003] Inter-chip interference (ICI) is a result of the multipath frequency 
selective channel in the CDMA downlink. The presence of ICI destroys the 

10 orthogonality of the Walsh spreading codes at mobile terminals. The 

challenge to the receiver design is even greater for a MIMO system in CDMA 
downlink. The receiver has to combat both the ICI and the co-channel 
interference (CCI) to achieve reliable communication. Therefore, 
interference cancellation at the mobile stations is an effective means of 

15 improving the receiver performance and link capacity. 

[0004] For a Single-In, Single-Out (SISO) or Single-In, Multiple-Out 
(SIMO) system, there are two types of linear equalizers: non-adaptive and 
adaptive equalizers. Non-adaptive equalizer usually requires matrix inversion 
and is therefore computationally expensive [1-4], especially when the 

20 coherence time of the channel is short and the equalizer requires frequent 
update. Adaptive equalizer is less computationally involved since it does not 
require direct matrix inversion, but it is not as robust as the non-adaptive 
equalizers since its performance is sensitive to the choice of parameters such 
as step size, initialization, etc [5-7]. To overcome these difficulties, in [8] an 

25 FFT-based linear equalization algorithm was proposed for SISO/SIMO 

systems that provides a good tradeoff between complexity and performance. 
The algorithm approximates the non-adaptive LMMSE algorithm by exploiting 
the banded Toeplitz structure of the correlation matrix of received signal 
vector. Another attractive alternative is the Kalman filtering approach 

30 proposed in [9] where an interesting state-space model is established to 
facilitate the application of the recursive Kalman filter. This approach 
outperforms the LMMSE approach at higher complexity. 
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[0005] For a Ml MO system, it has also been shown in [9] that both 
conventional LMMSE algorithm and the Kalman filter algorithm can be 
extended to the MIMO situation. In these MIMO solutions, both the ICI and 
CCI are suppressed in the linear equalization phase of the receiver. 
5 However, these solutions require complicated signal processing whose 
computational complexities may be beyond practical limits of the current 
hardware implementation. An FFT-based low-complexity MIMO LMMSE-FFT 
algorithm was recently proposed in [10] to overcome these difficulties. On 
the other hand, attempts to combine the non-linear decision feedback 
10 interference cancellation with the LMMSE equalization are also found in the 
literature, for example, [1 1]. However, these algorithms performs the decision 
feedback directly at the received signal, and thus requires the impractical 
assumption that all the active Walsh codes are known at the mobile receiver 
in order to reconstruct the transmitted chip sequences. 
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SUMMARY 

[0006] A system according to embodiments of this invention provides a 
multiple transmit antenna, multiple receive antenna (MIMO) receiver design 
for the communication downlinks such as those used in CDMA technology. 
5 Two algorithms referred to as the MIMO LMMSE-FFT and MIMO LMMSE-SIC 
(Successive Interference Cancellation) algorithms, are described in detail. In 
embodiments of the invention, the interference cancellation step is achieved 
without the impractical assumption of the knowledge of all the active Walsh 
codes in the systems. 

10 [0007] These and other features, aspects, and advantages of 

embodiments of the present invention will become apparent with reference to 
the following description in conjunction with the accompanying drawings. It is 
to be understood, however, that the drawings are designed solely for the 
purposes of illustration and not as a definition of the limits of the invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 is illustrative of a Ml MO-CDMA system. 

Figure 2 shows transmit signaling at antenna m of the system of figure 1 . 

Figure 3 is a block diagram showing the MIMO-LMMSE algorithm. 

5 Figure 4 is a block diagram showing a LMMSE-SIC receiver. 

Figure 5 is an illustration of a mixed traffic 1X EV-DV system. 

Figure 6 is a illustration of the conditional mean estimator. 

Figure 7 is a graph showing simulation results for the MIMO LMMSE-SIC 
algorithm. 

10 Figure 8 is a graph showing simulation results for the MIMO LMMSE-SIC 
algorithm with 2 iterations. 
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DETAILED DESCRIPTION 

[0008] Embodiments of this invention provide a multiple transmit 
antenna, multiple receive antenna (MIMO) receiver design for the 
5 communication downlinks such as those used in CDMA technology. 

[0009] Multiple transmit, multiple receive antenna MIMO systems offer 
potential for realizing high spectral efficiency of a wireless communications 
system. The information theoretic studies in [1] and [2] establish that in an 
independent at-fading channel environment, the capacity of such an MIMO 

10 system increases linearly with the number of antennas. In [3], a practical 
MIMO configuration, a Bell Labs Layered Space-Time (BLAST) system, is 
deployed to realize this high spectral efficiency for a narrow-band TDM A 
system. A simpler form of BLAST, the vertical BLAST (V-BLAST) is later 
proposed in [4, 5] to simplify the transceiver BLAST design. In V-BLAST, 

15 independent parallel data sequences are transmitted and a successive co- 
channel interference CCI cancellation/data detection algorithm is used for 
efficient multi-symbol detection. Note that the successive detection algorithm 
is also similar to those proposed in the CDMA multi-user detection literature, 
for example [6] [7]. 

20 [001 0] In the downlink of a CDMA like system, most of the research up- 
to-date focus on the application of the advanced signal processing techniques 
to offset the performance degradation due to the loss of Walsh code 
orthogonality caused by the frequency selective multipathchannel. Channel 
equalization is a promising means of improving the receiver performance in a 

25 frequency selective CDMA downlink. Current research encompasses two 
favors of linear equalization, namely non-adaptive linear equalization [8] [9] 
[10] [11] and adaptive linear equalization [12] [13]. Non-adaptive linear 
equalizers usually assumes "piece-wise" stationarity of the channel, and 
designs the equalizer according to some optimization criteria such as Linear 

30 Minimum Mean Square Error (LMMSE) or zero-forcing, which in general leads 
to solving a system of linear equations by matrix inversion. This can be 
computationally expensive, especially when the coherence time of the 
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channel is short and the equalizers have to be updated frequently. On the 
other hand, adaptive algorithms solve the similar LMMSE or zero-forcing 
optimization problems by means of stochastic gradient algorithms and avoid 
direct matrix inversion. Although computationally more manageable, the 
5 adaptive algorithms are less robust since their convergence behavior and 
performance depend on the choices of parameters such as step size. To 
overcome these difficulties, in [14] an FFT-based linear equalization algorithm 
was proposed that provides a good trade off between complexity and 
performance. The algorithm approximates the non-adaptive LMMSE 
10 algorithm by exploiting the banded Toeplitz structure of the correlation matrix 
of received signal vector. Another attractive alternative is the Kalman filtering 
approach proposed in [15], where an interesting 2 state-space model is 
established to facilitate the application of the recursive Kalman filter. This 
approach outperforms the LMMSE approach at a slightly higher complexity. 

15 [001 1] Applying MIMO configuration to the CDMA downlink presents 
significant challenge to the receiver design, as the receiver has to combat 
both the inter-chip interference (ICI) and the CCI in order to achieve reliable 
communication. It has been shown in [15] that both conventional LMMSE 
algorithm and the Kalman filter algorithm can be extended to the MIMO 

20 system. The attempts to combine the non-linear decision feedback 

interference cancellation with the LMMSE equalization are also found in the 
literature, for example, [16]. However, these algorithms performs the decision 
feedback directly at the received signal, and thus requires the impractical 
assumption that all the active Walsh codes are known at the mobile receiver 

25 in order to reconstruct the transmitted chip sequences. 

[0012] In embodiments of the invention, we propose a MIMO LMMSE- 
SIC algorithm combines the idea of linear equalization and non-linear 
decision feedback successive interference cancellation to improve the 
performance of a conventional linear MIMO LMMSE equalizer. One key 
30 advantage of the proposed MIMO LMMSE-SIC algorithm is that the 

interference cancellation step is achieved without the impractical assumption 
of the knowledge of all the active Walsh codes in the systems, unlike many 
other interference cancellation based algorithms found in the literature. 
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Noteworthy, in the proposed MIMO LMMSE-SIC algorithm, we have also 
incorporated the so-called conditional mean estimator to provide soft 
decisions in the decision feedback process. The simulation results suggest 
that the soft decisions alleviate the error propagation problem in the SIC 
5 associated with hard decision feedbacks. 

MIMO Signal Model for CDMA Downlink 

[0013] Consider an M transmit antenna, N receive antenna MIMO 
CDMA system as illustrated in Figure 1. Note that since Nokia's MIMO 
proposal on the transmitter side is yet to be finalized, in the receiver algorithm 
10 development process we have assumed a rather simple "serial to parallel 
split" transmit multiplexing, in order to make our receiver solutions general 
enough for all possible MIMO transmit multiplexing methods. The modulated 
symbol streams are split at the transmitter into M sub-streams before 
transmitted across the M transmit antennas. 

15 [0014] As shown in Figure 2, the signal model at the m t h transmit 

antenna is given as follows, assuming K active Walsh codes in the system: 

K 

dm{i) = c(i) ^ ^2 aekak, m (J) s k{i ~ jO) (1) 
fc=l m 

where m and /care chip, symbol, transmit antenna and spreading code 
indices. The base station scrambling code is denoted by c{i). Meanwhile, era 
20 stands for the power assigned to spreading code k (same for all antennas), 
ak.m is the information symbol sequence for spreading code k at antenna m 
and Sk is the km spreading code. Note that in this model we have implicitly 
assumed that the same set of Walsh codes are used across all the transmit 
antennas. 

25 [001 5] The transmitted signals propagates through the MIMO multipath 
fading channel denoted by Ho, ...,Hl, where each matrix is of dimension 
A/Ax/W, where A denotes number of samples per chip. The signal model at the 
receive antennas are thus given by the following equation, after stacking up 
the received samples across all the receive antennas for the /f/ichip interval. 
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y 4 = SMidi-i + m. (2) 

1=0 

Note that - y * UHi> .y*,iV"J j s Q f length A/ A, and each small vector y, ;n 
includes all the temporal samples within the ith chip interval. Meanwhile, L is 
the channel memory length, d *-* = (» ~ 0> -i<to(* - Q] T is the 
5 transmitted chip vector at time /' - /, and m is the (A/A) x 1 dimensional white 
Gaussian noise vector with n * ™ -A^CO, cr 2 I) Note tnat ^denotes noise 
variance and I is the identity matrix. Furthermore, in order to facilitate the 
discussion on the LMMSE receiver, we stack up a block of IF + 1 received 
vectors: 

Ji+Fx-F = H^i+F.i-F-L + n i+ jp :i _jp (3) 

10 

where 2F + 1 is the length of the LMMSE equalizing filter and 

y i+F:i - F = [yf +F ,...,yf_. J r] T , «2F + 1)JV A x 1) 

*k + F:i-F = K T +i r,...,Il|l F f J ((2F + 1)NA X 1) 

= [<% + f>-,<%-f-l\ T , ((2f + i + l)Afxl) 



?io • • • Hi 



((2F + 1) NA x (2F + L + 1)M) 



where the dimensions of the matrices are given next to them. Note that to 
keep the notation more intuitive, we keep the subscripts at a "block" level. For 
instance, yH-F-.i-F \ s the vector that contains blocks yi+Fj Ji-F where 
each block is a vector of size A/A * 1 . 



MIMO LMMSE Chip-Level Equalization 
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[0016] The block diagram of the MIMO receiver with chip-level 
equalizer is shown in Figure 3. It is essentially a straightforward extension of 
the LMMSE chip equalization initially designed for a SISO system [8]-[15]. 
After the chip-level equalizer, the orthogonality of the Walsh code is partially 
re-installed and all the desired symbols from each transmit antenna are 
detected with a simple code correlator which correlates to the desired 
spreading code. Note the descrambling process is also included in the code 

correlator. Defining an error vector of Z ~~ ^ Yi+f^-f anc j an error 

covariance matrix Rzz = £ t zz ^, the MIMO LMMSE chip-level equalizer W 
is the solution of the following problem: 

W^ £ - argtmn Traced) = arg min E\\d, - W*y i+M || 2 , (4) 

whose optimal solution is given by: 

= £7rfBL _1 Hi:i. (5) 

Where R = E ly i+F-.i-FY^-Fn-F] j S the correlation matrix of the 

2 

received signal, d is the transmitted chip power. Meanwhile, although 
H is fixed for a given channel realization and is not a function of symbol 
index /', here we use the notation H /+F:/+1 , H ;: , and H M: ,_ F L to represent 
the sub-matrices of H that are associated with d /+ F:/+i, d, and d i-vj-f-l in 
the expansion of the matrix-vector product: 



di 

_ di-i^-jp-i 

= H<+F:t+ld» + F:<+l + H«di 4- Hi_l : i_JP-£,di-l:<_p_l,. (6) 



Constrained Optimization Based Non-linear 
LMMSESIC 
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[0017] Both MIMO LMMSE and MIMO LMMSE-FFT belong to the 
category of linear equalization methods. In this section, we introduce the non- 
linear decision feedback to the MIMO LMMSE equalizer to improve the overall 
receiver performance. The resulting LMMSE-SIC (LMMSE- Successive 
5 Interference Cancellation) receiver is illustrated for a 2 Tx, N Rx MIMO 
system in Figure 4. There are two symbol-detection paths in the block 
diagram where the first/second path aims at detecting symbols transmitted on 
the first/second Tx antenna. The First path is similar to that of a conventional 
MIMO LMMSE receiver, where the chip-equalizer — using frequency domain 
10 representation for ease of exposition— is designed such that both the 

inter-chip interference (ICI) and the spatial co-channel interference (CCI) from 
the other transmit antenna is suppressed in the estimated chip sequence 

di{f), in order to facilitate the de-spreading and symbol detection by 
restoring the orthogonality between the Walsh codes. Note that only the 
15 symbols carried on the desired users' Walsh codes are detected. As shown in 
Figure 4, these detected symbols from the first path are directly fed-back and 
used in the symbol detection of the second path. 

[0018] The second path of the LMMSE-SIC receiver deviates from that 
of a conventional MIMO LMMSE receiver in that the chip equalizer W2(f) does 
20 not attempt to directly generate the chip sequence from the second transmit 
antenna. Instead, it attempts to generate a weighted sum of the chip 

sequences from both transmit antennas while suppressing all the ICIs, 

i a (f) = Mf) + Kid x {f) + Mf)- (7) 

i.e. 

this filter provides complete temporal ICI suppression while keeping a 
25 "controlled" residue CCI whose strength is denoted by a design 

parameter b 2 ,i. Furthermore, this design parameter is jointly optimized 
with the filter coefficients to achieve the best performance. The 
detailed derivation of this joint optimization is provided in the next 
section. Since the same set of Walsh codes are used in both Tx 
30 antennas, it follows that a simple despreading process gives us the symbol 
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estimate °" s ^ (for a particular Walsh code) that also includes a controlled 
residue CCI: 

&.(/) = «(/) + fciiaiC/) (8) 

where °"a ~~ ^^ a ^i anc ' ^ is the spreading length. It becomes obvious 

5 that assuming correct symbol estimations from the first path, the 

symbol estimate for the second Tx antenna is given by: 

Soft channel-coded bits can be easily extracted from the symbol estimates 

%(/)) «2(/) before being sent to the decoder to complete the receiver 
10 processing. As shown in a later section, significant performance gain can be 
achieved with LMMSE-SIC, compared to a conventional LMMSE receiver. 
Furthermore, the LMMSE-SIC algorithm achieve the performance gains 
without having to make the impractical assumption of knowing all the Walsh 
codes. This key advantage makes LMMSE-SIC especially attractive for a 
15 fixed voice-data system such as 1X EV-DV or HSDPA and the like, where the 
desired user usually accounts for only 10-50 percent of the overall transmit 
power in the cell. We illustrate such a mixed voice-data 1X EV-DV system in 
Figure 5, where the data user of interest and another data user each 
assumes about 20 percent of the transmit power, while the rest of the 
20 transmit power is assigned to voice users in the cell besides some 

housekeeping overheads such as pilot, synchronization, etc. Also note that in 
this example system, data users have a spreading length of 32 while the 
voice users have a spreading length of 64. 

Joint Equalizer/Feedback Weights Optimization 

25 [0019] We mentioned in the 2 Tx Ml MO example that the filter 

coefficients and the feedback weight b 12 requires joint optimization. In this 
section, we formulate this joint optimization problem for a general Mix MIMO 
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system and show that the problem can be converted into a series of M M VDR 
problems that are easily solved with the Lagrange Multiplier algorithm. 

[0020] For ease of exposition, we return to the matrix-vector notation 
used in equations (3 -6). The joint equalizer/feedback weights optimization for 
a MTx, N Rx MIMO system can be formulated as the following LMMSE 
problem with a lower-triangular structural constraint on the feedback weights: 

W^B^^argmmTraceCR^) = argmm 25||B"d, - W"y i+M || 2 , 



(10) 



[0021] Note that in this case the error vector is defined as 
z — B d* — W y»+jP:,-jr> A d j rect application of Lagrange multipliers to 
10 (10) proves to be difficult. To exploit the structure of the problem, we re- 
organize the structural constraint in (10) into two constraints that are more 
compact[17]: 



&m >m = 1, m = 1, ■ ■ ■ , M and 
Zh = 0 



(11) 



where ^ m m {^)mjm and 
15 in which b m is the m t h column of B. Furthermore, Z is defined as 



(13) 



where each sub-matrix Zm is an M * M diagonal matrix denoting the index of 
the null elements in b/, that is, (Zm)k,k = 1 if the km element in bmis constrained 
to be zero, and (Zm)k,k = 0 if it is not. Similarly, we define another diagonal 
20 matrix N = diagfM, ...,Nm} that denotes the index to the non-zero elements in 
b. It is easy to see that Zm + Nm = I by definition, where I is the identity matrix. 
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With these definitions, we set up the MMSE problem to provide optimal W 
and B with a generalized feedback structure: 

w * P * }B P P t = Qi[g ^nE\\B Il d i -W Ir y i+F:i ^\\ 2 > s.t. 

i) m ,m = 1, m = 1, ■ • - j M and 

Zh = 0. (14) 

[0022] To proceed, note that the cost function J in (14) can be re- 
5 written as: 

J = EWW H y i+F: ^ F -B H <Uf 

M M 

= E ||wSH i: ,-bS|| 2 +i:^Vw m (15) 

m=l m=l 

where V = ***** + a ^ w ™ is the nwi columns of W. Note that **** is 
a submatrix of H that is obtained by excluding H/:/ in the matrix H. Meanwhile, 
in the derivation we have also substituted y i+fj-f from (3). We now use the 
10 identity I = Zm + Nmto break up H/:/and bHmas: 

Hi : » = Hi : iAf m 4- H* :<J 2 m , 

= bfiAk + bJX. (16) 

[0023] Substituting (16) into (15) and after some algebra, we get 

J = E ll<^ - b^ m || 2 + E w^V m w m (17) 



where we have defined ^ m — ^:s^mH^ + V since ^""^"contains the 

15 controllable non-zero elements in m , we minimize the first term in J by 
setting 

b£A/™ =w^H i;i ^ m . (18) 
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Meanwhile, we observe that the second sum in J is optimized by individually 
minimizing 

w£f = aigminw^V m w m! s.t. w^li m = 1; for m = 1. M (19) 

[0024] Note that hm is the mm column of the matrix He/, and the 
h™ = 1 

constraint m ^ n follows from (18) and the original constraint bm.m = 1 . 
We have now decoupled the complex optimization problem (10) into a series 
of simple MVDR optimization problems (19) to get the filter coefficients wm. 
The non-zero elements in the feedback weight vectors bm are then obtained 
from (18). A straightforward application of the Lagrange Multiplier method to 
(19) leads to the final solution to the joint optimization 
problem: 



opt 



{\>°*) H M™ - (w£)'**H* t A/ r m , (20) 

and the elements of resulting error correlation matrix is given by 

(R* 2 )m,n = (w^)"(V + U i:i Z m 2 n U^) W ° n pt - (21) 

Soft Decision Feedback 

[0025] In the 2 Tx, 2Rx example illustrated in section 4, the estimated 

symbol from the first path is used in (9) to obtain the symbol estimates 

for the second path. If these symbol estimates (-^from first path are 
generated by hard-decisions, then the performance degradation for the 
second path can be severe due to error propagations. 

[0026] The error propagation problem can be alleviated by using the 
soft decision feedback. There are many types of soft estimators in literature, 
we found that the so-called conditional mean based estimator to be both 
analytically pleasing and easy to implement. This conditional mean estimator 
(also known as MMSE estimator) is also used in similar decision feedback 
algorithms such as [16]. Consider the general signal model for an arbitrary 
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symbol at the output of the de-spreader: 

T kt m{j) = ca km (j) + 71, (22) 

where c is a multiplicative factor and k,m, j are Walsh code, transmit 
antenna and symbol indices, respectively. Also note that we assume n 

to Gaussian: n ~ a «)and both c and °™ are functions of the 

spreading gain G, transmit power a k , noise level a 2 , channel H and the 

filter \Nopt. Once c and n are obtained, the conditional mean 
estimate of ak.mu) ' s given by 

a k , m (j) ^ E[a k>m (j)\r k>m (j)}, (23) 



which has a closed-form solution if we invoke the Gaussian assumption on n. 
For example, when the symbols of interest are QPSK modulated with unit 
energy, the solution reduces to a simple hyperbolic tangent function: 

[0027] The reason that the conditional mean estimator works well can be 
intuitively explained by observing Figure 6, where the real part of the symbol 

estimate ^( fl M»0)) j s plotted as 10 a function of the real part of the initial 

noisy soft estimate ^ e ( r *."»C?)) The hyperbolic function essentially acts as a 
clipper, which provides near-hard decision output when the initial estimate is 
\Re(r km (j))\ 

reliable, i.e. when is large. On the contrary, if the initial 

estimate is less reliable, the output is reduced to alleviate the possible error 
propagation impact. 

Optimization of Detection Order and Further Iterations 

[0028] In the LMMSE-SIC algorithm discussed above, we have .xed the 
detection order to be (1, 2, ...,M) rather arbitrarily. However, the performance 
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of the successive detection can be improved by optimizing the detection order 
[4] [5] [18]. In [4], the detection order for a similar V-BLAST problem is chosen 
so that the worst SNR among M data streams is maximized. In particular, let 
denote an arbitrary ordering and let Q be the set of all possible orderings. The 
5 cardinality of the set is | CI \ = M\ and the following problem is solved to obtain 
the optimal ordering p 

M 

/3 = axgmax mm. *y m (oj) , (25) 

where I'm ft") denotes the effective SNR associated with the detection of the 
mth transmitted data sequence for a given detection order . The effective 
10 SNR can be equivalently defined either at chip-level or at symbol level. The 
chip-level SNR for the mth data sequence is defined as a function of the error 
covariance matrix Rzz. 

which reduces to the following simple form considering the results in (20) and 
15 (21): 



7m (co) = cr*li2V-V)hm. (27) 

[0029] Once the effective SNR is defined, the problem in (25) can be 
solved by a "localized optimization" procedure similar to the one first proposed 
in [4]. In a localized optimization procedure, at each stage of the detection 
20 process, the sequence that has the best SNR is chosen for detection. The 
detail of this local optimization is not the focus of our discussion here. 

[0030] The performance of the LMMSE-SIC algorithm can be improved 
if we allow more complexity and introduce further iterations after all the data 
sequences are detected using the successive algorithm described above. 
25 Without loss of generality, we assume that in the .rst iteration the detection 
order is cu = (1, 2, • ■ -,M) and the filter matrix W and feedback matrix B is 
obtained by solving the optimization problem (10). In the second iteration, 
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however, the design of the matrix pencil W,B should be changed to reflect the 
fact that at the beginning of the second iteration, we have already obtained an 
initial estimate of all the transmitted data sequences. To this end, we change 
the constraint in (10) and rewrite the optimization problem as: 



W opt B op * 



axgmin£;||B H d,-W jr y i+F; ,_ F || 2 3 



fat A 



1 



(28) 



[0031] It is observed that in (28), none of the elements in the feedback 
matrix B is constrained to be zero. Therefore, it allows maximum amount of 
decision feedback in the spatial dimension and thus improve the effective 
SNR for each data sequence. Note that since (28) can also be transformed 
into the form shown in (14) with the compact constraint representation, the 
solution in (20-21 ) applies to (28) as well. However, it is easy to see that in 
this case, 

Z x = ■ ■ ■ = Z M = O, (29) 
and as a result 

Vi = ■ • • = V M - V. (30) 
[0032] Consequently, the solution to (28) assumes a simpler form: 
w qp* _ V lh ™ 
(*T) H = «)" rt H«. (31) 



[0033] We further denote lm as the effective SNR for the mth data 
sequence in the second iteration. It is easy to show that: 



(32) 
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[0034] Note that unlike Vm, V is independent of the detection order. 

Consequently, both w™, bm and ,m are independent of the detection order, 
meaning that all these quantities can be computed at the same time, unlike 
the local optimization process required in the first iteration. The detection 

5 order in the second iteration is thus obtained by simply sorting all the ,f ™ in 

the descending order. Furthermore, denoting lm as the SNR for mm data 
sequence in the first iteration, one can easily show that 

< (33) 
for any 1 * — m — ^ . Note that this inequality holds regardless of the 
10 detection order used in the first and second iterations. This inequality also 

signifies the benefit of having a second iteration in the LMMSE-SIC algorithm, 
provided that the added complexity is within the confine of a practical 
implementation. 

FFT-based Complexity Reduction for a 2 Tx LMMSE-SIC 
15 Algorithm 

[0035] In [19], an LMMSE-FFT algorithm originally proposed for a 
SISO/SIMO system in [14] is extended to the MIMO system. It was shown in 
[19] that the complexity of the an MIMO LMMSE algorithm can be greatly 
reduced with the FFT-based approach that exploits the block Toepltiz 

20 structure of the received signal correlation matrix R, and approximate R 1 with 
S \ where S is the associated block circulant matrix. It will be desirable if a 
similar lowcomplexity approach exists for the MIMO LMMSE-SIC algorithm as 
well. Unfortunately, since the matrices to be inverted in the LMMSE-SIC 
algorithm, namely Vm,m = 1, ■ ■ ,M, are not block Toeplitz, the direct 

25 extension of the FFT-based approach to the LMMSE-SIC algorithm is not 
possible for a general M Tx, N Rx MIMO system. For the special case where 
the number of transmit antenna M = 2, the matrices Vi and V2 can be written 
as a rank-one or rank-two update of a block Toeplitz matrix, which allows a 
low-complexity inversion that is similar to the MIMO LMMSE-FFT approach. 
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We focus on the .rst iteration of LMMSE-SIC and assume that the detection 

order is » = (*.»>. 

One can show that in this case 

V 2 = V, (34) 
and they relate to the correlation matrix R by 

V x - R-o>ibf 

V 2 = R-^H^Hg^R-^hihf-^hahf. (35) 
[0036] Substituting (35) into (20) and using the identity that 



(36) 



A -1 c _ (A — cr|cc ff )~ 1 c 
c jH A~ 1 c — c A (A — cr^cc lf )- 1 c' 

we can write the optimal filters as a function of correlation matrix R: 

Wl " hfVr 1 !!! "hfR-ihi' 

V2 X h 2 (R-gfojhf)- 1 ^ 
W2 - hfV^-hfXR-a^hf)-!^" (37) 

[0037] It is clear in (37) that the FFT-based approach can be used 
15 directly to get wi by approximating R with its associated circulant matrix 

S: R" 1 w S -1 .. On the other hand, the application of the FFT approach to W2 

requires a little more work. Observing that ^ ~~ rank-one update of 

the block Toeplitz matrix R, we invoke the matrix inversion Lemma [20] and 
obtain: 

20 ^-^ hlh ^ R +i_ 0 5hfR-ih 1 wS + i-o3hfs-im- (38) 

[0038] Note that in (38) we have already used the approximation 
R 1 a* S 1 From ( 37 _ 38 ) we have shown that for a 2 Jx Ml MO system 
using LMMSE-SIC algorithm, the filters wi and W2 can be obtained using FFT- 
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based approach to compute the inversion of the associated circulant matrix S, 
instead of the direct inversion of the correlation R. As a result, the complexity 
of the MIMO LMMSE-SIC algorithm for a 2Tx system can be close to that of 
the MIMO LMMSE-FFT algorithm, provided that no further iterations are used 
5 in the LMMSE-SIC. 
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Simulation Results 

[0039] The proposed MIMO LMMSE-SIC algorithms is evaluated in a 
realistic simulation chain, and the simulation parameters are tabulated in 
Table I. 



Parameter Name 


Parameter Value 


' System 


CDMA 1X/EVDV 


Spreading Length 


32 


Channel Profile 


Vehicular A 


Mobile Speed 


50 km/h 


Filter Length 


32 


Number of Tx/Rx Antennas 


2/2 


Modulation Format 


QPSK 


Information Data Elate 


312 kbps 


Turbo Code Rate 


0.6771 


Geometry 


6 


Number of Walsh Codes Assigned to the user 


3 


Total number of Active Walsh Codes in the system 


25 



Table 1: Simulation parameters. 

5 

[0040] Note that since the project is ongoing, the 14 results presented 
here are not complete and more results will be included in a future revision of 
this document. Nevertheless, these results provide good insights about the 
effectiveness of the algorithm proposed in this report. The benefit of the non- 
10 linear decision feedback is demonstrated in Figure 7, where the MIMO 

LMMSE-SIC algorithm is compared with the conventional LMMSE algorithm. 
Note that the impact of detection order optimization is included in these 
results. There are three LMMSE-SIC curves shown in this figure, representing 
the hard-decision feedback method, the conditional mean estimator based 
15 soft decision feedback method and the ideal feedback case. It is observed 
that even with the hard decision feedback, the MIMO LMMSE-SIC algorithm 
outperform the conventional LMMSE by about 1 dB at high SNR. Another 0.5- 
1dB can be gained if we replace the hard decisions with the soft decisions 
generated by the conditional mean estimator. It is important to note that these 
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gains are achieved with the assumption that the mobile receiver has only the 
knowledge of 3 out of 25 active Walsh codes. 

[0041] In Figure 8, the impact of further iterations in the LMMSE-SIC 
algorithm is investigated using the same example. Note that for both one 
5 iteration an two iterations case, the conditional mean estimator based soft 
decision feedback are used. In this case, adding a second iteration actually 
degrades the performance at high SNR. This detrimental effect, also known 
as "ping-pong" effect in the iterative detectors, is mainly attributed to the 
accumulation of decision errors in the iterative process. Obviously, for the 
10 MIMO LMMSE-SIC algorithm, further research is needed to reduce the 
impairment caused by the accumulative error propagation. 

[0042] Embodiments of the invention describe a MIMO receiver design 
that provides the best performance and complexity trade off. An embodiment 
of the invention is a potential candidate reference receiver for MIMO 
15 standards for 3GPP2. 

[0043] Although described in the context of particular embodiments, it 
will be apparent to those skilled in the art that a number of modifications and 
various changes to these teachings may occur. Thus, while the invention has 
been particularly shown and described with respect to one or more preferred 
20 embodiments thereof, it will be understood by those skilled in the art that 
certain modifications or changes, in form and shape, may be made therein 
without departing from the scope and spirit of the invention as set forth above. 
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ABSTRACT 

[0044] A system according to embodiments of this invention provide a 
multiple transmit antenna, multiple receive antenna (MIMO) receiver design 
for the communication downlinks such as those used in CDMA technology. 
Two algorithms referred to as the MIMO LMMSE-FFT and MIMO LMMSE-SIC 
(Successive Interference Cancellation) algorithms, are described in detail. In 
embodiments of the invention, the interference cancellation step is achieved 
without the impractical assumption of the knowledge of all the active Walsh 
codes in the systems, unlike many other interference cancellation based 
algorithms found in the literature. 
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